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We apply a modified mean-field density functional theory to determine the phase behavior of bi- 
nary mixtures of Stockmayer fluids whose spherical constituents interact according to Lennard- Jones 
(LJ) pair potentials with embedded pointlike dipole moments. On the basis of systematic numer- 
ical calculations we construct the global phase diagrams of these systems in the three-dimensional 
thermodynamic space of temperature, pressure, and chemical potential difference of the two com- 
ponents. The vapor-liquid, isotropic liquid - isotropic liquid, isotropic liquid - ferromagnetic liquid, 
and ferromagnetic liquid - ferromagnetic liquid first-order phase separations are investigated. The 
loci of the second-order isotropic fluid - ferromagnetic fluid transition are calculated from Landau 
theory. Liquid-vapor and liquid-liquid critical lines, tricritical lines, triple lines, and lines of critical 
end points of the binary Stockmayer mixtures are also determined. We discuss how the topology of 
the phase diagrams change upon varying the strengths of the two dipole moments of the two species 
as well as their sizes. 

PACS numbers: 64.10.+h, 64.75.+g, 68.18.Jk, 64.70.-p 



I. INTRODUCTION 

There are various classes of fluids whose structural and thermodynamical properties are strongly influenced by 
the presence of either magnetic or electric dipolar interactions between their constituents: polar molecular liquids, 
dipolar liquid crystals, ferrofluids, and magneto- or electrorheological fluids. These systems respond to external 
electric or magnetic fields and some may even exhibit long-ranged spontaneous orientational order in the liquid 
phases. The strongly anisotropic and long-ranged character of the dipolar forces causes particular difficulties for the 
theoretical description of these systems. By using simulations, integral equation theories, density functional theories, 
and phenomenological approaches significant progress has been achieved in understanding one-component dipolar 
fluids (see, e.g., Refs. Q and and references therein) albeit several aspects remain still under debate 

If isotropic phases of dipolar fluids form interfaces (such as liquid-vapor or fluid-wall interfaces) the broken trans- 
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lational invariance gives rise to interesting local orientational order tied to the interface [jl 
orientational profiles have become experimentally accessible through ellipsometry [9l. IToL [TlL 12t. L_ 
the orientational order at (critical) liquid-liquid 0, 0] and at (noncritical) fluid- vapor |ll , Il2( interfaces of several 
binary mixtures of polar (A) and nonpolar (B) components have been studied. In particular the growing spatial exten- 
sion of orientational order normal to the interface and proportional to the bulk correlation length upon approaching 
the continuous demixing transition of A-rich and B-rich phases in the bulk has been analysed. These data have been 
discussed in terms of universal scaling functions. Various forms of their shape have been prop osed 0, llfj| , including 
the occurrence of a surface correlation length different from the bulk correlation length |lll . These surprising 
interpretations lack systematic theoretical analyses, which are particularly important because the ellipsometry data 
provide only indirect access to orientational profiles so that theoretical guidance is required. It is our goal to provide 
such systematic theoretical studies for the aforementioned intcrfacial phenomena near continuous demixing and to 
^ ■ extend them to wetting phenomena. As a necessary first step the present work aims at providing an overview of the 
global phase behavior of polar binary liquid mixtures. This is an indispensable prerequisite for studying the intcrfacial 
properties in a second step. Moreover, in view of the numerous interesting features even of bulk properties of dipolar 
fluids HH3.Q- the analysis of dipolar mixtures is also important in its own right. 

Already some time ago there was an interesting first step towards the understanding of interfacial properties of 
dipolar liquid mixtures. In a combined experimental and theoretical study the liquid- vapor surface tensions of two 
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polar binary liquid mixtures have been studied for various temperatures and as function of composition |l3j . Using 
density functional theory (DFT) also relative adsorption as well as concentration and orientational profiles have been 
calculated [13 . By resorting to an improved DFT it is our intention to expand the above analysis towards focusing on 
liquid-liquid demixing transitions, the formation of phases with long-ranged orientational order, critical phenomena, 
and wetting phenomena and thus reaching out to the aforementioned more recent experimental studies. 

There is a sizeable body of theoretical research available which is devoted to the stud y of bulk properties of polar 
liquid mixtures. This encompasses Monte Carlo simulations Il4. lla . jla . ll7L fl8L [ill l20t l2ll 123 . l23j | . molecular dynamics 
simulations [TII2I I2H IH l2Sl. I29L integral equations IMEllIIIIl], cluster aggregation models |H HI IHUll , 
construction of equations of state |39j|. thermodynamic perturbation theory |40l l4ll lilj. structure formation in 
external fields 01 > an d van der Waals and Ising fluids . There are few theoretical studies concerning the structure 
formation of dipolar mixtures at interfaces such as DFT and Monte Carlo studies at charged walls [46ll47| and integral 
equatio n ap proaches for dipolar mixtures in random confined geometries |48l |49| . Beyond the aforementioned study 
in Ref. [L| the structure of dipolar mixtures at interfaces between fluid phases has been investigated by Monte Carlo 
simulations of the liquid-liquid benzene-water interface |5(Tj and by molecular dynamics simulations [Hlj and lattice- 
gas models |52j of the liquid- vapor interface of water-methanol mixtures. In the context of this body of research the 
present study is focused on obtaining global phase diagrams using an appropriate version of DFT which opens the 
possibility for interfacial studies later on. 

Models for dipolar fluids have to capture three basic features of the pair potentials between the fluid particles: 
short-ranged repulsion (which in general is anisotropic reflecting molecular shapes) , long-ranged attractive dispersion 
forces (containing an omnipresent isotropic contribution ~ r~ 6 at large distances r between the particles plus higher 
order anisotropic terms), and the dipolar interaction ~ r -3 . 

As long as one is not aiming for a quantitative description of a specific system but for general phenomena and trends 
the so-called Stockmayer model has turned out to be rather useful (see Refs. pj and y| and references therein). It 
considers spherical particles interacting with Lcnnard- Jones (LJ) potentials plus pointlikc permanent dipoles at the 
center. Such binary systems are characterized by three thermodynamic parameters (temperature T and two chemical 
potentials \ia and fis for the two species A and B), three LJ interaction energy parameters e.y (i, j = A,B), three 
LJ molecular diameter parameters <7y , and two dipole moments tua and tub ■ After choosing an energy scale for 
ksT = and the among the and a length scale for the number densities among the tr^ one is left with 
a six-dimensional interaction parameter space. By using DFT we analyze the three-dimensional (T,p, A/i) phase 
diagrams (where p is the pressure of the system and Afi = fig — ^a) as they change along certain paths within 
suitable subspaces of this interaction parameter space. 

This Stockmayer model provides an approximate model for dipolar molecular fluids as well as an effective model 
for ferrofluids with the solvent degrees of freedom taken into account via effective interaction energies e,j . Ferrofluids 
are often polydisperse [2ll l53l | . In this respect the study of binary Stockmayer models provides also insight into the 
effects of polydispersity for ferrofluids in the special case of bimodal distributions of their sizes or of their magnetic 
moments. 

Whereas one-component fluids with simple, radially symmetric pair potentials give rise to only two fluid phases, i.e., 
liquid and vapor, binary fluids exhibit typically three fluid phases: vapor, A-rich liquid, and B-rich liquid. The latter 
two emerge via phase segregation from a homogeneous liquid upon lowering the temperature. There is a line of critical 
points (with high temperatures) limiting a two-dimensional sheet of first-order liquid-vapor phase transitions, a line 
of critical points (at intermediate temperatures) limiting a two-dimensional sheet of first-order demixing transitions 
into ^4-rich and B-rich liquid, and a triple line of phase coexistence between all three fluid phases formed as the 
intersection of the two aforementioned sheets; at high temperatures the triple line ends at a critical end point where 
the line of critical demixing points meets the sheet of liquid- vapor coexistence (see, e.g., Rcf. |54j and references 
therein). At low temperatures freezing sets in such that there the triple line ends in a tetra point where the vapor 
phase, ^4-rich liquid, i?-rich liquid, and the solid phase coexist. (For the topology of this phase diagram see Ref. |55|^ 
In the present study we focus on sufficiently high temperatures so that freezing is not of concern. (For freezing in 
one-component dipolar fluids see Ref. 0] and references therein. F or q uasicrystalline order and glasslike structure 
formation in two-dimensional binary dipolar systems see Refs. |56ll57ll5q .) The introduction of dipole moments shifts 
these phase boundaries and may even give rise to the formation of new liquid phases with spontaneous long-ranged 
orientational order. 

In a recent series of papers jS^, l6fl l6ll | Range and Klapp have also extended the modified mean-field density 
functional theory to study the phase behavior of binary dipolar mixtures. There they have focused on the investigation 
of the isotropic fluid - ferromagnetic fluid phase separations. Studying dipolar hard sphere mixtures (and briefly 
mentioning binary Stockmayer fluid mixtures in Ref. |60j| ') they have found that the isotropic fluid - ferromagnetic 
fluid phase transitions are shifted towards higher densities as compared with the corresponding one-component dipolar 
fluids. These general trends are in agreement with our findings for the corresponding phase transitions in binary 
Stockmayer fluid mixtures. However, there are important differences to the case of dipolar hard sphere systems, 
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because due to the additional presence of isotropic dispersion forces Stockmayer fluids exhibit also a vapor phase, 
which is indispensable for studying the interfacial properties mentioned at the beginning of the Introduction. These 
dispersion forces modify also the other phase transitions. Therefore in the following we explore the wider phase 
behavior of binary Stockmayer fluid mixtures. In Sec. |il|we formulate the DFT for Stockmayer fluid mixtures. Our 
results are presented in Sec. Ill and summarized and discussed in Sec. IV. 



II. DENSITY FUNCTIONAL THEORY FOR STOCKMAYER FLUID MIXTURES 



A. Microscopic model 



We study binary Stockmayer fluids which consist of spherically symmetric particles interacting via Lennard- Jones 
(LJ) potentials with parameters aij and = A, B): 



wf/{r 12 ) = 4e, 
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In the following for the unlike LJ size parameter we assume cr^ 
interaction due to point dipolcs embedded at the centers of spheres of diameter an : 



(1) 



(an + crjj)/2. In addition there is the dipolar 



miirij 



[mi(wi) • r 12 ][m 2 (uj 2 ) • t"is 



[mi(o)i) ■ m 2 (w 2 )] 



0(ri2 - 



(2) 



where particle 1 (2) of type i (j) is located at r± (r 2 ) and carries a dipole moment of strength rrii (rrij) with an 
orientation given by the unit vector mi(wi) (m 2 (cj 2 )) with polar angles w\ = (9\, <f>\) (u> 2 = ($2, 4>2))\ r i2 = v\ — r% is 
the difference vector between the centers of the particle i and the particle j with 7*12 = |f 12I and O(r) is the Hcaviside 
step function. Thus the total pair potential is given by 



(3) 



In order to apply density functional theory, following Barker and Henderson |62l | this interaction potential is decom- 
posed into a short-ranged repulsive reference part 



w T ij (ria) = 9(f« - ri 2 )w ?;j (ri 2 ,cji,u) 2 ) = @(aij - r 12 )w\- } (r 12 ) 
and into a long-ranged attractive excess part 

Wij C (r 12 ,u)i,U2) = ©(V12 ~ (r^w^w-i). 



(4) 
(5) 



This decomposition allows one to choose a suitable reference system (a hard sphere system with temperature dependent 
diameters) for which reliable approximations for the free energy are known. The excess part w^ c (r\ 2l lui,lu 2 ) is treated 
perturbatively in a suitable way (see below). 



B. Modified mean-field density functional theory 



Our analysis is based on the following grand canonical variational functional, which is the binary extension of the 
one component functional used in Refs. and |6^|: 

y = fHs(PA,PB,T) + Pi [ ( W ) M 47ra i H] + + PlfX " ^ 

i=A,B %=A,B 

where V, T, and pi denote the volume of the fluid, its temperature, and the chemical potential of the ith component, 
respectively, f3 — l/(fcsT) with the Boltzmann constant Ub, Pi = N%/V is the number density of the species i taken 
to be spatially constant for the present bulk study. oti(u)) is the orientational probability density distribution (i.e., 
Pi(ui) = piOiiiuj) with J du>ai(u>) = 1 is the number density of particles of type i with orientation ui). The first term 
in Eq. JHJ is the free energy density of the reference system due to Mansoori et al. [6^ : 

r H i( PA , PB ,T) = /3- 1 J2 PiHAf^-iH/r 1 - 

i — ' 7T 

i=A,B 



I -Co 
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1-6 6(i-6) 2 
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where A4 denotes the thermal dc Broglie wavelength of the particles i and 

£fc = ~ E A- = 0,1, 2, 3 , (8) 



6 i=A,B 



with temperature dependent hard sphere diameters as given by Barker and Henderson |62l |: 

d u (T) = H' dr 12 (l - exp(-/3< e/ (r 12 ))). (9) 



For &AA = a BB Eq. (J7J reduces to the Carnahan-Starling formula. In Eq. © the second contribution to the grand 
canonical functional takes into account the loss of entropy if the orientational distributions of the particles are not 
isotropic. For the interaction contribution due to the long-ranged part of the pair potential between the spherical 
particles with embedded point dipoles we adopt 



a- 



E P*Pi J d3ri J J duJl j ^2ai(wi)a J (a;2)exp(-/3u'[/ / (r 12 ))/f/ c (ri 2 ,wi,a; 2 ) (10) 

i,j=A,B 

where 

/f/ c (ri 2 ,a;i,w 2 ) = 1 - cxp(-/3</ c (Vi 2 , u x , oj 2 )) (11) 

are the Mayer functions with the excess interaction potential. Equation I|1U|) amounts to using the low-density 
approximation for the pair distribution functions: 

^(r 12 )=exp(-/3 W [f(r 12 )). (12) 

For uniaxial molecules the orientational distribution functions oti(oS) depend only on the angle 9 and thus allow for 
the following expansion into Legendre polynomials: 



where 

2Z + 1 f 1 



27rai(w) = a, (cos (9) = ^ a\ l) Pi(cos6) (13) 



i° = / dxai(x)Pt(x). (14) 



Following Refs. jE3,E3] m order to calculate the spatial integrals in Eq. (|10J) we consider a finite volume V which has 
the shape of a prolate rotational ellipsoid with aspect ratio k (k > 1). The Mayer functions can be expressed in terms 
of the rotational invariants §1^1 (for the definition of these functions see, e.g., Ref. |(34j) so that the corresponding 
integrals yield 



i,j=A,B 1=0 



where 



^ g^. ^r^i roc p 

Wy = — -8i >1 I{k)m i m j - /T 1 ^ - ^ 3/2 J dr l2 r\ 2 J duj l du 2 f^ c {r l2 ,ui,u 2 ,ui 2 )^* m {u l ,u) 2 ,uj l2 ). (16) 

Further details of this expansion can be found in Ref. |64|. The function I{k) describes the shape dependence of the 
grand canonical potential, which enters only via the coefficients u[j and is given as 

^=s^-w^w^ k+v ^y <i7) 

Therefore the grand canonical functional reads 

y = f? e f(PA,pB,T) + k B T E Pi dxa i (x)ln[2a l (x)}+ ^ ^ u^af 'aj 5 - E ftMi- (18) 
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The equilibrium configuration for given T, /j,a, and \xb is determined by minimizing the total grand canonical functional 
with respect to pa, Pb-, an d the functions ola{x) and as(x): 



1 dfl 
V dp~ 



0. 



A,B 



(19) 



i sn 



V 6ati(x) 

Equations {ISJl and ([TSJ> lead to 



= 0, i = A,B. (20) 

a jlli {x),p i 



/I oo 
da;a i (aj)lii[2a i (s)]+2 ^ PjE^'W"' 1 ^ ' i = A,B. (21) 

1 j=A,B 1=0 



Equation lt2U|) yields 



ai(x) = Q exp ( -P^(2l + l)Pi(x) Pi^af ] , i = A,B, (22) 

{=1 3=A,B 



with 



1= /'dxcxpf + 1)/^) J] P^WM 



For the second term in Eq. 1|18[) one finds after some elementary calculations: 



i = A,B. (23) 



/l oo 
dxa,(x)ln[2« 4 (x)] = fc B Tln(2C l )-2 £ £ u^afaf , i = A,B. (24) 

Equations i|14[l and (|22|l lead to a system of coupled nonlinear equations for the coefficients aSJ± and aSg : 
(0 _ 21 + 1 /-i ^(*) CX P + WW E j= a,b Pi^S a f\ 

f\ dzcxp (-^E,= l(2i + l)i^(«) Ej=A,B Pj^S^ 



Equation (|21|l together with Eqs. (|23|l and (|24() give the following expressions for the chemical potentials of the two 
components: 

IM = p? s + k B Tln(2Q) + ^ Pi u i?> i = A,B. (26) 

J=A,B 

The chemical potentials pa and ps can be eliminated from Eq. I|18|) and therefore fl can be expressed as 
o °° 

~=f%0>A,PB,T)- J2 Pitf S (PA,PB,T)- J2 PiP^M^f- ( 27 ) 

i=A,B i,j=A,B 1=0 

Due to Q = —pV the corresponding pressure of an equilibrium phase as a function of the densities and the temperature 
is 

oo 

p = -f r H e s f ( PA ,PB,T)+ Pl p? s (PA,p B ,T)+ PiPjJ2 u< S a i )a f- ( 28 ) 

i=A,B i,j=A,B 1=0 

These results show that the knowledge of the functions u^(T) is necessary to obtain the coefficients oq of the 
oricntational distribution functions (Eq.{5SJ)) and to obtain the different thermodynamic functions (Eqs. (^26128^ . 
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In computing these functions we use the same approximations as in Ref. [64| . i.e., we expand the Mayer functions 
in E q. (Iltjf) for small dipole moments which allows us to perform the angular integrations analytically. As in 
Ref. [§4j we have determined the expansion coefficients up to 0(m 8 ). This gives a satisfactory accuracy for reduced 
dipole moments m* = rrii/iJ caao"aa ^ 2 |64| . The omission of the higher order terms implies a truncation in the 
summations at I = 4 in all equations which contain such a summation over I. The corresponding coefficients ufj(T) 
are: 

,(o) 



"u 



^ / dri 2 rl 2 (l - exp(-^w^ J (r 12 )) - ^^m?m|7 4 (,9ey) - ^g-mfmpw^eij), (29) 

P JO •> rr , 7 - ;,rr ,, 

ug> = -ym imj /(fc) - l|^ m 3 m 3 /7(/?£y . )5 (30) 



(2) _ 2 o , _ 327T/33 __ 4 __ 4 

" 375^ "'<"V ;; 6125a, 



T m l m j h(Pe lj ) - 9 m.mjho^eij), (31) 



(3) = 167T/3- ^ 
y 25752(7,;, 



T mimp 7 {ptij), (32) 



and 



with 



Uy = ~ 99225a 9 - TOi;Wj 10 ( ^ 



dxx~ n exp[4y(s 



a:-"-!-")]. (34) 



C. Phase equilibria and phase diagrams 

The phase diagrams at a given temperature T follow from requiring the equality of the chemical potentials of both 
components and of the pressure for the densities pp and pp characterizing the coexisting phases I and II: 

P\p[ I \& l l I \x) = P\p["\a[ II Hx)- 

The functions a^P (x) and af 1 ^ 1 (x) (i = A, B) denote the corresponding equilibrium orientational distributions ob- 
tained from Eqs. (|25|) and (|22|l . As mentioned in the Introduction we consider four kinds of two-phase equilibria: 
isotropic liquid - isotropic vapor, isotropic A-rich liquid - isotropic B-rich liquid, isotropic liquid - ferromagnetic 
liquid, and ferromagnetic A-rich liquid - ferromagnetic B-rich liquid. For binary fluids Gibbs' phase rule allows for 
three-phase coexistence lines and four-phase (tetra) coexistence points in the thermodynamic space (T,p, A/i). For 
isotropic phases the corresponding equations for the chemical potentials and the pressure are simpler because for them 
oti(x) = 1/2 (i.e., af ) = 1/2 and af = if I > 1). 



In order to locate the surface of critical points separating the isotropic and the ferromagnetic liquid phases by 

HIS 

(0 



s sep ; 

second-order phase transitions, similar as in Ref. Ml we expand the orientation dependent ideal gas contributions 

to the interaction part flint of the grand canonical functional for small deviations from isotropy, i.e., for small a. 
Using the results for the one component dipolar fluid with I = 1 the leading orientation dependent contribution to 
the grand canonical potential turns out as 



i,j=A,B V 
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where Sij is the Kronccker delta symbol. Obviously, for given (3, pa, and pb the isotropic configuration minimizes the 

grand canonical functional if the quadratic form defined by Eq. H^fiJ) is positive definite. Provided pau aa + 2/ (3/3) > 0, 
this quadratic form is positive definite if and only if its determinant is positive |66l | . For all the cases studied it turns 
out that this inequality is satisfied. Thus the zero of the determinant defines the surface which separates the isotropic 
liquid phase from the ferromagnetic liquid phase: 

(pAU^a + {^PBU { g B + - pAPB (u^b) = 0. (37) 

In the one-component limit Eq. (|37|l reduces to the equation for the critical line given in Ref. [64[ ■ 

III. RESULTS AND DISCUSSION 

In the following we shall use reduced quantities: T* = ksT/tAA as reduced temperature, p* = (pa + Pb)&aa 
as reduced total density, Xi = pi/{pA + Pb) as the concentration of species i, p* = pa AA /eAA as reduced pressure, 
p* = {pi — ksT ln(Af / a ' AA )) 1 6 aa as the reduced chemical potential of species i, and m* - nxij \J ^aa^aa as the 
reduced dipole moments. For our numerical analysis we consider the following choices for the Lennard- Jones potential 
parameters: cbb/caa = 0.95, cab/^aa = 0.75, and <jaa = &bb = gab- The case o aa ^ obb will be discussed in 
Subsec. Ill C. The bulk phase diagrams of the binary dipolar fluid mixtures are calculated in the thermodynamic 
space spanned by temperature T , pressure p, and the difference of the chemical potentials Ap = pb — PA of the two 
components. In order to obtain the phase boundaries in this space for fixed T the corresponding phase equilibrium 
equations (Eq. 1(35(0 have to be solved simultaneously under the constraint 

Ap-p B (p I A ,p I B ,T) + p A {p A ,PB,T)=0. (38) 

The pressure p follows from Eq. 

A. Binary mixtures of particles with equal size and equal dipole moment 

As a first step we have determined the phase diagrams of non-polar binary Lennard-Jones fluid mixtures (m A = 
m* B = 0) for the aforementioned parameters. Based on our systematic calculations in the thermodynamic space 
(T,p,Ap) these phase diagrams take on shapes as shown schematically in Fig. ^a). Si and 5*2 are sheets of first- 
order phase transitions separating the vapor phase from the fluid phase and the A-rich liquid phase from the B-rich 
liquid phase, respectively. These sheets are bounded by lines L\ and L 2 , respectively, of second-order phase transitions 
marking the onset of liquid- vapor separation (L\) and demixing (L 2 ) upon lowering the temperature. The intersection 
of Si and 52 forms the triple line TL of three-phase coexistence ending at the critical end point CEP where L 2 hits 
Si . This type of phase diagrams corresponds to type II in the classification scheme of phase diagrams of binary fluid 
mixtures given by Scott and van Konynenburg [671 l68j|. For the non-polar binary fluid mixture Fig. [21 presents our 
numerical results for the pressure-temperature (Fig. Ufa)) and the density-temperature (Fig. Efb)) phase diagrams for 
Ap* = 0.3. At low temperatures (T < T tr ) the A-rich liquid coexists with the vapor. At the triple point temperature 
T t * — 1.003 three phases coexist: vapor, B-rich liquid, and A-rich liquid. Above T tr up to the liquid-vapor critical 
temperature (T c = 1.153) there is coexistence between the vapor and the B-rich liquid and between the A-rich 
liquid and the B-rich liquid. In the temperature range X'j' 1 '' 1 < T < Tc with J^"* 1 * = 1.429 there is only coexistence 
between the A-rich liquid and the B-rich liquid. 

Figure^b) shows the schematic phase diagram of the binary liquid mixture with the same LJ parameters but with 
m A = m* B = 1 . The occurrence of a ferromagnetic phase gives rise to the appearance of the sheet S3 which separates 
the isotropic liquid at high temperatures from the ferromagnetic liquid at low temperatures. The nature of this sheet 
S3 differs from that of Si and S2 in that S3 is the locus of second-order phase transitions. In the thermodynamic 
space (T,p,Ap) the surface S3 is given by the simultaneous solution of Eqs. (|37|) and (|38(l . (With the knowledge 

of pa, Pb, and T Eq. 1)28(1 gives the corresponding pressure p, using the fact that on S3 af^ = 1/2 and = 
for I > 1.) The intersection of S3 and Si is a line Mi of critical end points. Along this line the ferromagnetically 
critical liquid phase coexists with the isotropic vapor phase. The other line M 2 of critical end points is given by the 
intersection of the sheets S2 and S3. The sheet S3 divides the triple line into two parts. The upper part describes 
isotropic vapor - isotropic B-rich liquid - isotropic A-rich liquid three-phase coexistence, while the lower part describes 
an isotropic vapor - ferromagnetic B-rich liquid - ferromagnetic A-rich liquid coexistence. Both critical lines Li and 
L 2 are in the isotropic part of the thermodynamic space (T,p, Ap). Figure [3] shows the corresponding results of our 
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FIG. 1: Schematic bulk phase diagrams of binary Stockmayer fluid mixtures in the thermodynamic space of temperature T, 
pressure p, and chemical potential difference Afi — fi B — of the two components for various equal dipole moments. The 
phase diagram in (a) represents a non-polar (m A — m B = 0) binary fluid mixture which exhibits three phases: vapor, A-rich 
liquid, and B-rich liquid. They are separated by sheets Si and S2 of first-order phase transitions, delimited by lines L\ and 
Li, respectively, of second-order phase transitions. Si and S2 intersect along the triple line TL which ends at the critical end 
point CEP. The phase diagram in (b) represents a dipolar (m% = m B = 1) binary fluid mixture for which a third sheet 
S3 (hatched) appears, which separates the isotropic liquid from the ferromagnetic liquid. S3 is a sheet of second-order phase 
transitions. The intersection of the sheets Si and S3 forms a line Mi of critical end points. The intersection of the sheets S2 
and S3 forms another line M2 of critical end points. Upon increasing the dipole moments (m A = m B = 1.5) in (c) the sheet S3 
has been raised above L\ and has merged with Si into the sheet Si U S3 with a tricritical line (TC-line), which separates the 
first-order transition (lower part of the sheet) from the second-order transition (hatched upper part). The isotropic liquid- vapor 
transitions are reduced to two winglike surfaces, which are connected to the sheet Si U S3 along the triple lines TL\ and TL[ . 
These triple lines end in critical end points CEP\ and CEP[. The liquid- vapor critical lines L\ and L' ± are the upper edges of 
the wings. Upon a further increase of the dipole moments (rn* A = m B = 2) one obtains the phase diagram shown in (d) where 
there are only two sheets: S2 and Si U S3. The wings associated with the isotropic liquid - isotropic vapor transitions have 
disappeared and only the isotropic vapor - ferromagnetic liquid and the ferromagnetic A-rich liquid - ferromagnetic B-rich 
liquid phase transitions remain. 



numerical calculations for a binary Stockmayer fluid mixture with dipole moments m* A = m* B = 1 for A/i* = 0.3. The 
intersection of the plane A/j,* = 0.3 with the sheet £3 yields the critical (dotted) line in Fig. Ol^a). Its intersection with 
the first-order coexistence curve between vapor and ferromagnetic A-rich liquid is the critical end point CEP with 
T CEP = 1.082, p CEP = 0.03, and p CEP = 0.74 at Afj* = 0.3. Below Tcep the isotropic vapor phase coexists with 
the ferromagnetic A-rich liquid. The critical line of second-order phase transitions divides the A-rich liquid into an 
isotropic A-rich liquid phase and into a ferromagnetic A-rich liquid phase. In the temperature range Tcep < T < T tr 
with T t * = 1.12 the isotropic vapor coexists with an isotropic A-rich liquid phase. Above the triple point temperature 
there is also coexistence between the isotropic B-rich liquid and the isotropic A-rich liquid. This first-order phase 
transition ends in the critical temperature Tc* = 1.404. The critical temperature of the isotropic vapor - isotropic 
B-rich liquid coexistence is Tc h '^* = 1.296. 

Upon increasing the dipole moments the phase diagram shown in Fig. ^b) changes in two respects. First, the 
sheet S3 is raised due to the strengthening of orientational ordering. Second, close to the line Mi the character of the 
ferromagnetic - isotropic phase transition on S3 becomes first-order such that a line of tricritical points (TC) emerges 




FIG. 2: Phase diagrams of a non-polar {m\ = m* B = 0) binary Lennard- Jones fluid mixture in the pressure-temperature (a) 
and density-temperature (b) plane for A^i* = 0.3; ess/eAA = 0.95, cab/caa = 0.75, oaa = &bb = tab- T tr , Tc lv \ and T c 
denote the triple, liquid-vapor critical, and liquid-liquid critical temperature, respectively. In (a) all lines are lines of first-order 
phase transitions. The vapor phase occurs only at very low pressure. In (b) the full curves enclose two-phase regions. The 
dashed line indicates three-phase coexistence. Below the triple point temperature the vapor coexists with the A-rich liquid. 
Between Ttr and tJ there are three possible phases: the vapor, the isotropic B-rich liquid, and the A-rich liquid. Between 
Tc^ and Tc there is only B-rich liquid - A-rich liquid coexistence. At high densities the system freezes. The corresponding 
solid phases are not shown because they are not accessible by the present theory. 




FIG. 3: Phase diagrams of a binary Stockmayer fluid mixture for m\ = m* B = 1 in the pressure-temperature (a) and density- 
temperature (b) plane for A/U* = 0.3. The LJ parameters are the same as in Fig. |5] The dotted line denotes a line of 
second-order phase transitions. It ends at a critical end point (CEP) where the ferromagnetically critical liquid coexists with 
the non-critical vapor. The critical line divides the A-rich liquid region into an isotropic and ferromagnetic part. For more 
details see the caption of Fig. [5] 



from Mi and separates the sheet S3 into a continuous part (back part at high p and hatched) and a first-order part 
(front part at lower p and unhatchcd). This implies that M\ turns into a triple line TL*. Upon further increase of 
rriA and m# this sheet S3, divided by the TC- line, raises finally above the line £2- Moreover the upper part of Si, 
i.e., above TL* disintegrates into a left and right part leaving a gap between them. This means that TL* is split into 
TLi and TL[ and Li into Li and L[ in Fig. ^c). This also generates two critical end points CEPi and CEP[ which 
are raised above CEP. Figure Wic), corresponding to m* A = m* B = 1.5, displays this latter structure with Si U S3 
denoting the sheet which emerged from the lower part of Si in Fig. ^b) and S3 whereas the winglike structures in 
Fig. (He) are the remnants of the upper part of Si in Fig. IHb). They are connected to the sheet Si U S3 along the 
triple lines TLi and TL[ and bounded by the liquid- vapor critical lines L\ and L[. The isotropic vapor - isotropic 
fluid - ferromagnetic liquid triple lines TL\ and TL\ end in critical end points CEPi and CEP[. S2 and Si U S3 
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FIG. 4: Projections onto a plane p — const of the triple lines and of the liquid-vapor critical lines of a binary Stockmayer fluid 
mixtures with m* A — m* B = 1.5; the LJ parameters are the same as in Figs. |2]and|21 In accordance with Fig. Q^c) along TL 
there is isotropic vapor - A-rich ferromagnetic liquid - B-rich ferromagnetic liquid three phase coexistence. This triple line ends 
at a critical end point (CEP). Both TL\ and TL[ describe isotropic vapor - isotropic fluid - ferromagnetic fluid three-phase 
coexistence. These two triple lines end at critical end points CEP\ and CEP[. The dashed lines L\ and L[ represent the 
projections of the liquid- vapor critical lines delimiting the upper end of the wings shown in Fig. 0]c). Note that as in Fig. 
Au increases towards the left. 



T 1 



2.2 












(a) 


isotropic fluid 


(b) 




2 




v 


isotropic fluid 


trc .-' — 




T 1 " 




T (lv) 


/ \ ferromagnetic - 
/ \ liquid 


1.8 




ferromagnetic liquid 


vapor 




1.6 










1.4 




Au*= 4.0 


■ i A/- 4.0 

i.i.i 









0.5 1 1.5 2 
P 


0.2 0.4 


0.6 0.8 



1.8 j 



- 1.6 



FIG. 5: Phase diagrams of a binary Stockmayer fluid mixture for m* A = m* B = 1.5 (compare Fig. Wic)) in the pressure- 
temperature (a) and density-temperature (b) plane for A/i* = 4; the LJ parameters are the same as in Figs. EE] Below 
the triple point temperature Tt r the isotropic vapor coexists with the ferromagnetic liquid. Between Tt r and Tc lv ^ there are 
three possible phases: the isotropic vapor, the isotropic liquid, and the ferromagnetic liquid. The first-order phase transition 
between the isotropic fluid and ferromagnetic liquid turns into a second-order phase transition at the tricritical temperature 
Ttrc- Above the tricritical temperature the dotted line represents the corresponding line of critical points. For the value of the 
chemical potential difference A/x* = 4 there is no A-rich liquid - B-rich liquid phase separation (i.e., the plane A/x* = 4 does 
not intersect the sheet S2; see Figs. 0and^J. For more details see the caption of Fig. [5] 



intersect along the triple line TL which ends at the critical end point CEP. 

Figure 01 shows the corresponding results for the projections of the triple lines and of the liquid- vapor critical lines 
of a binary Stockmayer fluid mixture with dipole moments m* A = m* B = 1.5 onto a plane p = const. The triple line 
TL describes the isotropic vapor - ferromagnetic A-rich liquid - ferromagnetic B-rich liquid three-phase coexistence 
with the critical end point temperature T^ EP = 1.614. The triple lines TL\ and TL\ associated with the two winglikc 
surfaces (see Fig. ^c)) describe the isotropic vapor - isotropic fluid - ferromagnetic A-rich liquid and isotropic 
vapor - isotropic fluid - ferromagnetic B-rich liquid three-phase coexistences with the critical end point temperatures 
T* Ep , = 1.764 and T* EPi = 1.747. 

Figure [S] presents our numerical results for a binary Stockmayer fluid mixture with dipole moments m* A = m* B = 1.5 
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FIG. 6: Phase diagrams of a binary Stockmayer fluid mixture with m\ = m B = 1.5 (compare Fig. 0|c)) in the pressure- 
temperature (a) and density-temperature (b) plane for A/x* = 0.4; the LJ parameters are the same as in Figs. _______ For this 

choice of A/z, below the isotropic vapor - ferromagnetic A-rich liquid - ferromagnetic B-rich liquid triple point temperature Ttr 
the isotropic vapor coexists with the ferromagnetic A-rich liquid. Between Ttr and there are three possible phases: the 
isotropic vapor, the ferromagnetic B-rich liquid, and the ferromagnetic A-rich liquid. In (b) the two-phase region of A-rich and 

B-rich liquid coexistence ending at T c is rather narrow. Between and the tricritical temperature Ttrc the isotropic vapor 
coexists with the B-rich ferromagnetic liquid. Above the tricritical temperature Ttrc the transitions between the isotropic and 
ferromagnetic fluid are continuous and their loci are given by the dotted critical line. The dash-dotted line in (b) indicates the 
thermodynamically unstable isotropic liquid - vapor coexistence. 
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FIG. 7: Phase diagrams of a binary Stockmayer fluid mixture with m* A = m* B — 2 (compare Fig. [__[d)) in the pressure- 
temperature (a) and density-temperature (b) plane for Afi* = 0.5; the LJ parameters are the same as in Figs. |__E_| Below 
the triple temperature T tr the isotropic vapor coexists with a ferromagnetic A-rich liquid. The vapor phase occurs at such 
small values of p and p that it is not visible on the scales shown in (a) and (b). Within the very narrow temperature range 
Tr <T < T { c ll) there are three possible phases: the isotropic vapor, the ferromagnetic B-rich liquid, and the ferromagnetic 
A-rich liquid. Above the ferromagnetic B-rich liquid - ferromagnetic A-rich liquid critical temperature there is only 

coexistence between the isotropic fluid and the B-rich ferromagnetic liquid. Above the tricritical temperature Ttrc these phase 
transition are continuous forming the dotted critical line. 



for A/i* = 4. For this choice of A/x at the triple point T tr the following three phases coexist: vapor, isotropic fluid, 
and ferromagnetic liquid. At low temperatures T < T tr with T t * = 1.768 the ferromagnetic liquid coexists with 
the isotropic vapor. Above the triple point temperature, two isotropic fluids coexist at low and moderate densities, 
becoming identical above the liquid-vapor critical temperature yj'")* — 1.788. (The corresponding critical pressure 

and density are p£ = 0.185 and pc* = 0.286, respectively.) At higher densities, the isotropic liquid and magnetic 
liquid are separated by a first-order phase transition which becomes second order above the tricritical temperature 
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T* rc = 1.954. (The corresponding tricritical pressure and density are p\ rc = 0.964 and p\ rc — 0.579, respectively.) 
The critical line of second-order phase transitions divides the liquid range into an isotropic liquid phase and into a 
ferromagnetic liquid phase. At this value of the chemical potential difference there is no A-rich liquid - B-rich liquid 
demixing phase separation. 

Figure shows our numerical results for a binary Stockmayer fluid mixture with dipolc moments m A = m* B = 1.5 
for A/i* = 0.4. At the triple point temperature T± r the coexisting phases are: vapor, ferromagnetic B-rich liquid, and 
ferromagnetic A-rich liquid. At low temperatures T < T tr with T t * = 1.512 the isotropic vapor phase coexists with the 
ferromagnetic A-rich liquid. Above the triple point temperature, an isotropic vapor coexists with a ferromagnetic B- 
rich liquid and a ferromagnetic B-rich liquid coexists with a ferromagnetic A-rich liquid. The two ferromagnetic liquid 
phases become identical beyond the liquid- liquid critical point 7c with Tc = 1.694. (The corresponding critical 
pressure and density are pi — 16.83 and p£ — 1.17, respectively, at A/i* = 0.4.) Above the ferromagnetic A-rich 
liquid - ferromagnetic B-rich liquid critical temperature Tc the isotropic liquid and ferromagnetic B-rich liquid are 
separated by a first-order transition which turns second order above the tricritical temperature T t * c = 1.918. (The 
corresponding tricritical pressure and density are p\ rc = 0.994 and p* trc = 0.568, respectively, at A/i* = 0.4.) At 
temperatures T > T trc one finds second-order ferromagnetic phase transitions as given by the critical line. 

Fig. ^d) shows that by further increasing the dipole moments {m* A = m B = 2) the winglike surfaces associated 
with the isotropic liquid - isotropic vapor transition disappear and only the two sheets S2 and Si U S3 remain. 
Accordingly there is only one triple line TL, which describes the isotropic vapor, ferromagnetic A-rich liquid, and 
ferromagnetic B-rich liquid three-phase coexistence and which ends at the critical end point CEP where Li hits 
S\ U S3. Since the sheet S2 bounded by L2 lies beneath the sheet Si US3 all first-order liquid- liquid phase separations 
involve ferromagnetic phases. 

Figure \7\ shows our numerical results for a binary Stockmayer fluid mixture with dipole moments m A = m B = 2 
for A/i* = 0.5. The topology of Fig. 0is the same as in Fig. El only the numerical values for the triple, critical, 
and tricritical points are different. Below the triple point temperature T* r = 1.961 (with p^ r = 0.002 at A/i* = 0.5) 
an isotropic vapor coexists with a ferromagnetic A-rich liquid. The ferromagnetic liquid-liquid critical point is given 
by Tc = 1-98, Pc* — 7.59, and p c — 1.2, while the tricritical point is given by T t * c = 3.17. p* trc = 2.06, and 
p\ rc = 0.528. 

B. Binary liquid mixtures of non-polar and polar particles of equal size 

In this subsection we present our results for phase diagrams of non-polar [m* A = 0) - polar (m B 7^ 0) binary 
mixtures. The corresponding schematic phase diagrams are shown in Fig. |H1 In order to provide a convenient visual 
comparison Fig. l%Ta) shows again as a reference case the schematic phase diagram of a non-polar binary Lennard- Jones 
fluid mixture (m A = m* B = 0) with the parameter set ^bb/^aa = 0.95, cab/^aa = 0.75 and <jaa = <Jbb = &ab- 
Figure |S{b) shows the schematic phase diagram of the binary liquid mixture with the same LJ parameters but with 
m A = and m* B = 1. Similar to Fig. ^b) the occurrence of a ferromagnetic phase gives rise to the appearance 
of the sheet S3 of second-order phase transitions which separates the isotropic liquid phase from the ferromagnetic 
liquid phase. The ferromagnetic phase appears only in the B-rich liquid phase and at low temperatures while the 
A-rich liquid phase remains isotropic even at low temperatures. This behavior reflects the polar (B) and non-polar 
(^4) character of molecules. In accordance with this asymmetry the sheet S3 is tilted, ends at S2 and thus does not 
cut through the sheet S2. This leads to the line M2 of critical end points. The intersection of S3 and Si forms another 
line Mi of critical end points. The triple line (TL) is given by the intersection of Si and S2. The line Mi of critical 
end points divides the triple line into two parts. The upper part of the triple line describes isotropic vapor - isotropic 
B-rich liquid - isotropic A-rich liquid three-phase coexistence as in the case of symmetric dipolar mixtures with dipole 
moments m A = m* B = 1 (see Fig. H^b)). The lower part of the triple line describes isotropic vapor - ferromagnetic 
B-rich liquid - isotropic A-rich liquid coexistence. 

Figure [5] shows our numerical findings for a binary Stockmayer fluid mixture with dipole moments m A = and 
rn* B = 1 for A/i* = —0.27. Figure Eta) demonstrates that for this value of the chemical potential difference the plane 
A/i* = —0.27 cuts through the upper part of the triple line TL in Fig. OJb) so that the corresponding triple point 
T tr describes an isotropic vapor - A-rich isotropic liquid - B-rich isotropic liquid coexistence with T t * = 1.16 and 
p\ r = 0.101. The intersection of this plane with the sheet S3 yields the dotted critical line. The intersection of the 
plane A/i* = -0.27 with the line Mi in Fig. ED» gives the critical end point CEP where an isotropic vapor coexists 
with a ferromagnetically critical B-rich liquid. The thermodynamic parameters of this critical end point CEP are 
given by T£ EP — 1.07, p* CEP = 0.043, and p* CEP = 0.74, respectively, for A/i* = —0.27. Below the temperature 
Tcep the isotropic vapor phase coexists with the ferromagnetic B-rich liquid phase. Between the temperatures Tqep 
and T tr two phases coexist: the isotropic vapor and the B-rich isotropic liquid. Above T tr up to the liquid-vapor 
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FIG. 8: Schematic bulk phase diagrams of binary Stockmayer fluid mixtures in the thermodynamic space of temperature T, 
pressure p, and chemical potential difference A/i = p B — pa of the non-polar (m A — 0) and polar (m B ^0) components. For 
reference purposes the phase diagram in (a) represents the same non-polar (m A — m B — 0) binary fluid mixture as in FigHTal. 
The phase diagram in (b) represents a non-polar - dipolar (m* A = 0, m* B — 1) binary fluid mixture for which in comparison 
with (a) a third sheet S3 (hatched) appears, which separates the B-rich isotropic liquid from the B-rich ferromagnetic liquid. 
53 is a sheet of second-order phase transitions. Mi and M2 represent lines of critical end points at the intersection of S3 with 
Si and S2, respectively. In (a)-(d) the A-rich liquid is always isotropic. Upon increasing the dipole moment of the component 
B (m A = 0, m B = 1.5) in (c) the tilted sheet S3 has moved above L\ and L2 and has merged with Si and S2 into the sheet 
S1US2US3 with a tricritical line (rC-line), which separates first-order transitions (lower part of the sheet) from second-order 
transitions (hatched upper part). The loci of isotropic liquid - vapor transitions are reduced to two winglike surfaces, which are 
connected to the surface S1US2US3 along the triple lines TL\ and TL[. These triple lines end in critical end points CEP\ and 
CEP[. The liquid- vapor critical lines L\ and L[ are the upper edges of the wings. Upon a further increase of the dipole moment 
of component B (m A = 0, m B = 2) one obtains the phase diagram shown in (d) where there is only one sheet S1US2US3 and 
one wing attached to it. The surface S1US2US3 denotes the isotropic vapor - ferromagnetic liquid transition, while the wing 
is associated with the A-rich isotropic liquid - isotropic vapor transition. 



critical temperature (T c = 1.19) coexistence is possible between the vapor and the A-rich isotropic liquid phase and 
between the A-rich isotropic liquid and the B-rich isotropic liquid. In the temperature range 

T (lv) <T< T (U) w . th 

yW* _ 23 coexistence is possible only between the isotropic A-rich and B-rich liquid phases. (The corresponding 
critical densities are pi * = 0.31 and p c l = 0.55, respectively.) 

Figure I1UI shows our numerical results for a binary Stockmayer fluid mixture with dipole moments m\ = and 
rn* B = 1 for A/i* = —0.54. Figure ITUT a) displays that for this value of the chemical potential difference the plane 
A/i* = —0.54 cuts through the lower part of the triple line TL in Fig. OJb) so that the corresponding triple point T tr 
describes an isotropic vapor - A-rich isotropic liquid - B-rich ferromagnetic liquid coexistence with T t * = 1.05 and 
Pt r = 0.045. Below the triple temperature T tr the vapor phase coexists with the B-rich ferromagnetic liquid phase. 

Above T tr up to the liquid-vapor critical temperature (Tc = 1.21) coexistence is possible between the vapor and 
the A-rich isotropic liquid phase and between the A-rich isotropic liquid and the B-rich ferromagnetic liquid. The 

critical pressure and density of liquid-vapor coexistence are p c lv ^* = 0.122 and p[ } = 0.285, respectively. The A-rich 
isotropic liquid - B-rich ferromagnetic liquid phase coexistence extends to higher temperatures. The corresponding 
liquid-liquid critical point is not attainable at physically relevant fluid densities. (At high densities the system freezes, 
but the corresponding solid phases are not captured within the framework of the present theory.) 
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P P 

FIG. 9: Phase diagrams of a non-polar - polar (m\ = 0, m B = 1) binary Stockmayer fluid mixture (compare Fig. Elb)) in 
the pressure-temperature (a) and density-temperature (b) plane for Afj,* = —0.27. The LJ parameters are the same as in Fig. 
[5] Below the critical end point (CEP) temperature Tqep the isotropic vapor coexists with the B-rich ferromagnetic liquid. 
Between Tcep and Ttr the isotropic vapor coexists with the B-rich isotropic liquid. Between T tr and there are three 

possible phases: the vapor, the A-rich isotropic liquid, and the B-rich isotropic liquid. Between Tc and Tc only A-rich 
isotropic liquid - B-rich isotropic liquid coexistence is possible. The dotted line denotes a line of second-order ferromagnetic 
phase transitions. It ends at a critical end point CEP. 




FIG. 10: Phase diagrams of a binary Stockmayer fluid mixture for m* A — and m* B — 1 (compare Fig. IHIb)) in the pressure- 
temperature (a) and density-temperature (b) plane for A^i* = —0.54; the LJ parameters are the same as in Fig. [5] Below the 
triple point temperature Ttr the isotropic vapor coexists with the B-rich ferromagnetic liquid. Between Ttr and Tc'^ there are 
three possible phases: the vapor, the A-rich isotropic liquid, and the B-rich ferromagnetic liquid. The A-rich isotropic liquid - 
B-rich ferromagnetic liquid coexistence does not end in a critical point within the range of fluid densities. 



Upon increasing the dipole moment of component B the phase diagram shown in Fig. E[b) changes considerably. 
First, the sheet 5*3 moves up due to the strengthening of the dipolar interaction. Like in the case of symmetrical 
dipolar mixtures (m* A = m* B = 1.5, see Fig. ^c)) close to the line Mi the ferromagnetic-isotropic phase transition on 
the sheet S3 becomes first order. This means that a line TC of tricritical points emerges from Mi and separates the 
sheet 5*3 into a continuous part and a first-order part. In turn this implies that Mi turns into a triple line TL*. In the 
next step due to the rise of the sheet S3 (divided by the TC-line) the line M2 rises, too, and merges with the line L2 
so that S3 and S2 form a single sheet S2US3. Upon a further increase of ms the upper part of Si disintegrates into a 
left and a right part leaving a gap between them, where Si merges with S2US3 forming the surface S1US2US3. The 
remnants of the upper part of the sheet Si appear as asymmetrically winglike structures, where the triple lines TLi 
and TL'i stem from the triple line TL* which has split while the split line Li generates the liquid-vapor critical lines 
L\ and L[ (see Fig. [Etc)). The intersections of the corresponding triple and critical lines give two critical end points 



15 




FIG. 11: Phase diagrams of a binary Stockmayer fluid mixture for m* A = and m B = 1.5 (compare Fig. |SJc)) in the 
pressure-temperature (a) and density-temperature (b) plane for A/t* — 1; the LJ parameters are the same as in Fig. Below 
the tricritical temperature T trc the B-rich isotropic fluid coexists with the B-rich ferromagnetic fluid. This first-order phase 
transition turns into a second-order phase transition at the tricritical temperature T± r c- Above T trc the dotted line represents 
the line of second-order phase transitions. 




FIG. 12: Phase diagrams of a binary Stockmayer fluid mixture for m* A — and m* B — 1.5 (compare Fig. EIc)) in the pressure- 
temperature (a) and density-temperature (b) plane for A/i* = —2; the LJ parameters are the same as in Fig. At this value 
of the chemical potential difference only the A-rich isotropic liquid can coexist along a line of first-order phase transitions with 
the B-rich ferromagnetic liquid. The liquid-liquid coexistence does not end in a critical point within the range of physically 
relevant fluid densities. 



CEP\ and CEP[. FigurcOUc) corresponds to m* A = and m* B = 1.5 and displays this final structure with S1US2US3 
denoting the surface which has emerged from the sheets Si, S2, and S3 (shown in Fig. Hl^b)) upon the increase of the 
dipole moment m^. 

Figure 1111 shows our numerical results for a binary Stockmayer fluid mixture with dipole moments m* A = and 
m* B = 1.5 for A/i* = 1. For this value of A/i* the figure displays a vertical intersection of the phase diagram shown 
in Fig. |HJc) between the two critical end points CEP\ and CEP[ and close to CEP\. At low temperatures the 
B-rich isotropic fluid and the B-rich ferromagnetic fluid are separated by a first-order phase transition which becomes 
second order above the tricritical temperature T t * c = 2.03. (The tricritical pressure and density are p* rc — 1.82 and 
Ptrc = 0.65, respectively.) The dotted critical line corresponds to a cut through the second-order part of the surface 
SiUS 2 US 3 . 

Figure 1121 shows our numerical findings for a binary Stockmayer fluid mixture with dipole moments vn* A = and 
m* B = 1.5 for A//* = —2. This figure displays a vertical cut through the phase diagram in Fig. QJc) between CEP\ and 
CEP{ but close to the critical end point CEP[ . For this value of the chemical potential difference the A-rich isotropic 
fluid can coexist with the B-rich ferromagnetic fluid. Note that the B-rich isotropic fluid and the A-rich isotropic 
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FIG. 13: Phase diagrams of a binary Stockmayer fluid mixture for m* A = and m* B = 1.5 (compare Fig. [Etc)) in the pressure- 
temperature (a) and density-temperature (b) plane for A/j,* = —6; the LJ parameters are the same as in Fig. Below the 
triple point temperature Ttr the isotropic vapor coexists with the B-rich ferromagnetic liquid. Between T tr and there are 
three possible phases: the vapor, the A-rich isotropic liquid, and the B-rich ferromagnetic liquid. The A-rich isotropic liquid - 
B-rich ferromagnetic liquid coexistence does not end in a critical point within the range of physically relevant fluid densities. 

fluid are not separated by a demixing transition; in Fig. Hfc) both are in front of SiUS^US^. The liquid-liquid critical 
point lies outside the physically relevant fluid densities. 

Figure 1131 shows our numerical results for a binary Stockmayer fluid mixture with dipole moments m* A = and 
m B = 1.5 for A/i* = —6. In this case the plane A/i* = —6 cuts through the wing on the right side of the surface 
S1IIS2US3 in Fig. OJc). For this choice of An at the triple point T tr the following three phases coexist: vapor, A-rich 
isotropic liquid, and B-rich ferromagnetic liquid. At low temperatures T < T tr with T t * = 1.27 (p* r = 0.1) the B-rich 
ferromagnetic liquid coexists with the isotropic vapor. Above the triple point temperature at low densities the vapor 
coexists with the A-rich isotropic liquid, whereas at high densities the A-rich isotropic liquid coexists with the B-rich 
ferromagnetic liquid. The vapor and the A-rich isotropic liquid merge at the liquid-vapor critical point Tc Lv ^* = 1.34, 
p* = 0.13, and p* c = 0.28. The A-rich isotropic liquid - B-rich ferromagnetic liquid critical point occurs outside the 
physically relevant fluid densities. 

Figure[IJd) shows that by further increasing the dipole moment of the component B {m* A = 0, m* B = 2) the winglike 
structure on the left hand side in Fig. OUc), which is associated with the B-rich isotropic liquid - vapor transitions, 
disappears so that there remain only the two sheets SiUS^US^ and the wing on the right hand side associated with 
the A-rich isotropic liquid - vapor transitions. Accordingly there is only one triple line TL[, which describes the 
isotropic vapor - A-rich isotropic liquid - B-rich ferromagnetic liquid three-phase coexistence. It ends at the critical 
end point CEP[ where L[ hits the surface .SiUS^US^. Since the topology of the phase diagram in Fig. Etd) is the 
same as that of Fig. Efc) (only the left winglike surface in Fig. |Hfc) is missing), the topologies of different cuts at 
given chemical potentials are in agreement with the topologies of the phase diagrams in Figs. 1111131 Therefore we 
refrain from presenting our detailed numerical results for Stockmayer mixtures with dipole moments m\ = and 
m* = 2. 

C. Binary liquid mixtures of polar particles of different size 

In this subsection we present phase diagrams of binary Stockmayer liquid mixtures of polar particles (m A y^m B ) of 
different size. Since in the Stockmayer interaction potential the size parameters appear in the LJ contribution to the 
potential, we first study the effect of varying the ratio &b b / &aa on the phase diagrams of binary L J fluid mixtures 
(m A = m* B = 0) while keeping the energy parameters of the LJ potential fixed. 

Figure^lshows the projections of the triple lines of binary LJ fluid mixtures with different particle size ratios onto 
a plane p = const. These triple lines describe the vapor - A-rich isotropic liquid - B-rich isotropic liquid three-phase 
coexistence. The curves show that the sheet S2, which separates the A-rich and B-rich liquid phases (see Fig. DJa)) 
is increasingly tilted upon increasing the size ratio cfbb/o~aa- The critical end point temperature increases and the 
chemical potential difference Afi^, EP at the end point turns more negative, i.e., it is shifted towards the region of 
fluid rich in the small particles. This means that for large size ratios three phase coexistence occurs for a wider range 
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FIG. 14: Projections of the vapor - A-rich isotropic liquid - B-rich isotropic liquid triple line TL (see Fig. of binary 
Lennnard- Jones fluid mixtures {m* A = m* B = 0) onto a plane p = const for ebb/caa = 0.95, bab/caa = 0.75, and various size 
ratios obb I o aa- The triple lines end at different critical end points (CEPs, see Fig. 0. Note that as in Figs. Qand|l]A/i 
increases towards the left. 




FIG. 15: Phase diagrams of the binary Lennard- Jones fluid mixtures in the pressure-temperature (a) and density temperature 
(b) plane for different size ratios obbI^aa and Afi* . The LJ energy parameters are the same as in Fig. 1141 In (a) all liquid- 
vapor coexistence curves separate the vapor phase from the B-rich isotropic liquid phase. The liquid-liquid coexistence curves 
separate the B-rich isotropic liquid phase from the A-rich isotropic liquid phase. In (b) the full, dashed-dotted, and dotted 
curves enclose two-phase regions whereas the dashed lines indicate three-phase coexistence. The line code in (b) corresponds 
to that in (a) with the same choices of the size ratio obbIoaa and A[x* . 



of thermodynamic parameters. 

Figure ITSl represents our numerical results for cuts through the phase diagrams of binary LJ mixtures for different 
size ratios o~bb/o~aa an d chemical potential differences A/i*. For the presentation we have changed the choice for 
AyU* along with gbb/o'aa because there is no single value Ap* for which the triple points exist for the different values 
o~bb /o~aa=1, 1.05, and 2 1 / 3 (see Fig. 11411 . Figure ITST a) demonstrates that the liquid-vapor critical temperatures 
depend barely on the size ratio (T^* = 1.153, T 7 ^ = 1.153, and T C 3 * = 1.152) while the liquid-liquid critical 
temperatures exhibit a stronger size ratio dependence (T^f = 1.258, T^p* = 1.131, and T^* = 1.148). The 
corresponding critical densities are given by p^* = 0.273, p&* = 0.217, p^* = 0.195 and p^i = 0.708, p^ = 
0.502, PcP* = 0.453. The corresponding critical pressures of the liquid-vapor and A-rich liquid - B-rich liquid critical 

points arc p^* = 0.113, = 0.091, = 0.083 and = 1.418, = 0.219, p ( c f* = 0.139, respectively. 

Figurc lTSl shows that the A-rich liquid - B-rich liquid coexistence range within the pressure-temperature plane becomes 
narrower while within the density-temperature plane it grows wider with increasing size ratio. 

We now turn to binary Stockmayer fluid mixtures whose components have different dipole moments and particle 
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FIG. 16: Schematic bulk phase diagrams of binary Stockmayer fluid mixtures in the thermodynamic space of temperature 
T, pressure p, and chemical potential difference A/x = [ib — fJ-A for various dipole moments and size parameters. The phase 
diagram in (a) represents a binary dipolar fluid mixture with the parameters m A = 0.75, m B — 1.5, and <tbb/<taa = 1. The 
topology of this phase diagram is the same as that of a non-polar - dipolar binary fluid mixture with the parameters m* A — 0, 
rn% = 1.5, and obbIoaa = 1 (compare Fig. |5fc)). Upon increasing the size of the component B such that ubb/^aa = 2 1 / 3 
(m* A — 0.75, m B — 1.5) one obtains the phase diagram shown in (c) in which the sheet S3 appears again, which separates 
the B-rich isotropic liquid from the B-rich ferromagnetic liquid. The topology of this phase diagram is the same as that of 
non-polar - dipolar binary fluid mixture with the parameters m A — 0, m B = 1, and o BB /a AA = 1 (compare Fig. ffi b)). The 
phase diagram in (b) represents a binary dipolar fluid mixture with the parameters m A = 1, m* B — 2, and obbIoaa = 1. 
The topology of this phase diagram is similar to the topology of the phase diagram of non-polar - dipolar binary fluid mixture 
with the parameters m* A = 0, m* B = 2, and obbIoaa = 1 (compare Fig. |Ht c )), with the difference that at the bottom of the 
winglike structure attached to the surface S1US2US3 an A-rich ferromagnetic liquid phase appears. In this phase diagram a new 
sheet of second-order phase transition S3 appears, which separates the A-rich isotropic liquid from the A-rich ferromagnetic 
liquid. The intersection of the winglike surface and the sheet S3 forms a line M\ of critical end points. The intersection of the 
sheets S1US2US3 and S3 forms another line M2 of critical end points. Upon increasing the size of the component B so that 
cbb/ctaa = 2 1 / 3 {m* A = 1, m* B = 2) one obtains the phase diagram shown in (d) where a winglike surface of first-order phase 
transitions appears on the left part of the surface Si US2US3. This winglike surface separates the vapor phase from the B-rich 
isotropic liquid phase. 



sizes. This model gives a reasonable description of bidisperse ferromagnetic fluids and as such provides also a first 
approximation of actual polydisperse ferrocolloids. For these fluids the magnetic dipole moments of the particles 
are proportional to the particle volume, i.e., rrn oc af i7 i = A,B. In the following we restrict our study to binary 
Stockmayer fluid mixtures with dipole moment ratio m* B jm* A = 2, which implies the size ratio ubb jo AA = 2 1 / 3 . This 
motivates our choice of the size ratio used for the above study of the phase behavior of the binary LJ fluid mixtures. 

For the dipole moment ratio m* B / m* A = 2 in the following we compare the phase diagrams of binary Stockmayer 
fluid mixtures composed of particles with cither equal or different sizes. The corresponding schematic phase diagrams 
are shown in Fig. 1161 Figure ITBT a) displays the schematic phase diagrams of a binary Stockmayer fluid mixture with 
dipole moments m* A = 0.75, m B = 1.5, and with particle size ratio ubb/^aa — 1; the LJ energy parameters are the 
same as in Fig. El The topology of this phase diagram is the same as that of shown in Fig. IHtc). This fact is not 
surprising because in both cases m* B = 1.5 and the increase of m\ from m* A = to m A = 0.75 causes only a small 
change in the dipolar interaction energy relative to the LJ interaction energy. 
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FIG. 17: Phase diagrams of a binary Stockmayer fluid mixture for dipole moments m* A = 0.75 and m B = 1.5, and with size 
ratio obbIoaa = 1 in the pressure-temperature (a) and density-temperature (b) plane for Afi* = —6 (compare Fig. E2 a ))- 
The LJ energy parameters are the same as in Fig. 1141 Below the triple point temperature Ttr the isotropic vapor coexists 
with the B-rich ferromagnetic liquid. Between Ttr and Tc there are three phases: the vapor, the A-rich isotropic liquid, and 
the B-rich ferromagnetic liquid. The A-rich isotropic liquid - B-rich ferromagnetic liquid coexistence does not end in a critical 
point within the range of fluid densities. 
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FIG. 18: Phase diagrams of a binary Stockmayer fluid mixture for dipole moments m\ = 0.75 and m B = 1.5, and with size 
ratio obbIoaa = 2 1 ' 3 in the pressure-temperature (a) and density-temperature (b) plane for A/i* = —1.36 (compare Fig. 
HGf c) - ). The LJ energy parameters are the same as in Fig. 1141 Below the triple point temperature the vapor coexists with the 
A-rich isotropic liquid. Between Ttr and there are three phases: the vapor, the B-rich isotropic liquid, and the A-rich 
isotropic liquid. Between To and T c there is only vapor - B-rich isotropic liquid coexistence. 



Figure IT7I shows our numerical results for this binary Stockmayer fluid mixture with dipole moments m A = 0.75 
and m* B = 1.5, and with particle size ratio obbI&aa = 1 for A/i* = —6. Figure nTT a) demonstrates that the plane 
A/t* = —6 cuts through the wing attached to the right hand side of the surface S1US2US3 in Fig. I16f a). In accordance 
with the aforementioned topological similarity the phase diagrams shown in Fig. 1171 differ from those shown in Fig. 

IT31 only quantitatively. In Fig. El the liquid-vapor critical point is given by T^ lv) * = 1.384, pi lv) * = 0.138, and 
p(f v )* _ 0.277. The corresponding triple point temperature and pressure are T t * = 1.257 andpj r = 0.077, respectively. 

For the same dipole moments (m 1 ^ = 0.75 and m* B = 1.5) increasing the particle size ratio from <tbb/&aa = 1 to 
(Jbb/&aa = 2 1 / 3 the phase diagram shown in Fig. I16f a) changes considerably. The resulting phase diagram is shown 
in Fig. Iltjf c). The topology of this phase diagram is the same as in Fig. |Sfb). As discussed in Subsec. Ill B the phase 
diagram shown in Fig. |Hfb) emerges from the phase diagram shown in Fig. IHJc) by decreasing the dipole moment 
m* B . Figure iToT a') and Fig. HoT c) show that the same topological change in the phase diagrams can be attained by 
increasing the corresponding particle size ratio. Figure ITdT c'I tells that the dipole moment strength m* A = 0.75 is 
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FIG. 19: Phase diagrams of a binary Stockmayer fluid mixture with dipole moments m A = 1 and m* B — 2, and size ratio 
cbb/^aa — 1 in the pressure-temperature (a) and density-temperature (b) plane for A/i* = 5.5 (compare Fig. Iltjf b')'). The LJ 
energy parameters are the same as in Fig. 1141 Below the tricritical temperature Ttrc the B-rich isotropic fluid coexists with 
the B-rich ferromagnetic fluid. This first-order phase transition turns into a second-order phase transition at the tricritical 
temperature Ttrc- Above T trc the dotted line represents a line of second-order phase transitions. 




FIG. 20: Phase diagrams of binary Stockmayer fluid mixture for dipole moments m\ = 1 and m B = 2, and with size ratio 
obb/oaa = 2 1//3 in the pressure-temperature (a) and density-temperature (b) plane for Afi* = 5.5 (compare Fig. EJd)). 
The LJ energy parameters are the same as in Fig. 1141 Below the triple point temperature Ttr the isotropic vapor coexists 
with the B-rich ferromagnetic liquid. Between T tr and Tc lv ^ there are three phases: the isotropic vapor, the B-rich isotropic 
liquid, and the B-rich ferromagnetic liquid. The line of first-order phase transition between the B-rich isotropic fluid and B-rich 
ferromagnetic liquid turns into a line of second-order phase transition (dotted line) at the tricritical temperature Tt rc - 



insufficient to lead to the formation of an A-rich ferromagnetic liquid phase, so that for negative values of A/x* the 
phase remains isotropic, as in the case of non-polar - polar mixtures (see Fig. [ffib)). 

Figure ITHl shows our numerical findings for a binary Stockmayer fluid mixture with dipole moments m A = 0.75 
and m* B = 1.5, and with particle size ratio <jbb/&aa = 2 1//3 for A/z* = —1.36. Figure ITST a) demonstrates that 
the plane A/i* = —1.36 cuts through the upper part of the triple line TL in Fig. I16f c) so that the corresponding 
triple point T tr describes a vapor - A-rich isotropic liquid - B-rich isotropic liquid coexistence with T t * = 1.1 and 
p* tr = 0.04. The vapor-liquid and liquid-liquid critical points are given as Tc* = 1.258, p£ = 0.096, pc = 0.204 
and Tg = 1.18, pc = 0.126, pc = 0.462, respectively. The topology of these phase diagrams is the same as 
those of shown in Fig. [3 

Increasing the dipole moments from m\ = 0.75 and m* B = 1.5 to m* A = 1 and m* B = 2 at a fixed particle size ratio 
pbb/paa = 1 the phase diagram shown in Fig. I16f a1 changes considerably leading to the phase diagram displayed 
in Fig. I16r b'). This increase of the dipole moments causes the disappearance of the winglike structure attached to 
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the left part of the surface 1S1US2US3 (see Fig. H6f b)L which is associated with the B-rich isotropic liquid - vapor 
transitions. Moreover the winglike structure attached to the right part of the surface SiUS 2 US 3 also changes due to 
the appearance of the sheet $3 . This sheet separates the A-rich isotropic liquid at high temperatures from the A-rich 
ferromagnetic liquid at low temperatures. S' 3 is the locus of second-order phase transitions and is the remnant of the 
sheet ^3 (see Fig. d^)), which together with Si and £2 has formed the surface S1US2US3. 

Figure 1191 displays our numerical results for a binary Stockmayer fluid mixture with dipole moments ra\ = 1 , 
m* B = 2, and with particle size ratio obb/vaa = 1 for Ap* = 5.5. For this value of Ap* the figure shows a vertical 
slice through the left part of the surface SiUS^US^ in the phase diagram displayed in Fig. I16f b) . The topology 
of these phase diagrams is the same as those of shown in Fig. ^2 The tricritical point is given by T t * c = 3.28, 
p* trc = 2.538, and p* trc = 0.565. 

Increasing the particle size ratio from <jbb/&aa = 1 to cfbb/^aa = 2 1 / 3 for fixed dipole moments m* A = 1 and 
m* B = 2, the phase diagram shown in Fig. I16f b) transforms into that shown in Fig. I16f d). Upon the increase of 
the particle size ratio on the left part of the surface S1US2US3 in Fig. I16f d) a winglike structure appears, which is 
associated with B-rich isotropic liquid - vapor phase transitions. Through a similar topological change, the phase 
diagram shown in Fig. Efd) transforms into the phase diagram shown in Fig. Etc) upon a decrease of the dipole 
moment m* B from m* B = 1.5 to m* B = 1. 

Figure 1201 displays our numerical results for a binary Stockmayer fluid mixture with dipole moments m* A = 1 , 
rn* B = 2, and with a particle size ratio o~bb/o~aa = 2 for Ap* = 5.5. For this choice of Ap* Fig. 120( a) represents 
a vertical slice through the left part of the surface S1US2US3 in the phase diagram shown in Fig. I16f d~) where the 
attached winglike structure appears again. The corresponding liquid- vapor critical point, the tricritical point, and 
the triple point are given by T* = 1.71, p* = 0.09, p* = 0.142, T*. c = 1.77, p* trc = 0.37, p* trc = 0.3, and T*. = 1.64, 
p* r = 0.065, respectively. The topology of the phase diagrams shown in Fig. [201 is the same as of those of shown in 
Fig. m 

IV. SUMMARY 

In the present study of the phase behavior of binary Stockmayer (Eqs. (d~©) fluid mixtures the following main 
results have been obtained: 

(1) We have extended the one-component version of the modified mean-field density functional theory to that 
of binary Stockmayer fluid mixtures of particles A and B with different dipole moments, sizes, and Lennard- Jones 
interaction energies. In particular, we have focused on the explicit formulation of the grand canonical functional, the 
equilibrium pressure p, and the chemical potentials pa and pb of the two components [see Subsec. II B and Eqs. 

(2) This system exhibits 6 distinct fluid phases: isotropic vapor, isotropic fluid, A-rich isotropic liquid, A-rich 
ferromagnetic liquid, B-rich isotropic liquid, and B-rich ferromagnetic liquid. At a given temperature T the first-order 
phase transitions between them follow from the equalities of the chemical potentials and of the pressure between the 
coexisting phases [see Subsec. II C and Eq. I|35|l]. The loci of the second-order phase transitions between isotropic 
and ferromagnetic fluids can be determined analytically from a suitable Landau theory [see Subsec. II C and Eq. 

63] ■ 

(3) In the thermodynamic parameter space (T, p 1 Ap = pb — Pa) we have constructed the schematic phase diagrams 
of binary fluid mixtures of particles with equal size and equal dipole moment (see Fig. d- We have discussed the 
topological changes of these phase diagrams due to an increase of the particle dipole moments. The various shapes of 
the three-dimensional schematic phase diagrams are supported by numerically calculated cuts through them, which 
are two-dimensional phase diagrams at fixed chemical potential differences Ap. For dipole moments vtia = mB = 
Fig. El shows a cut through the phase coexistence surfaces (Fig. d a )) of first-order liquid- vapor and A-rich liquid - 
B-rich liquid phase transitions of Lennard- Jones fluid mixtures which are characterized by energy parameters saa, 
£bb, £ab and diameters <taa, cbb, &ab = {&aa + &bb)/2- For reduced dipole moments m* A = m B = 1 with 
in* = rrii/ eaa<J aa , i = A,B, Fig. EJdisplays a cut through the phase coexistence surfaces (Fig. db)) describing the 
vapor-liquid, the vapor - A-rich ferromagnetic liquid, and the A-rich isotropic liquid - B-rich isotropic liquid phase 
transitions of binary Stockmayer fluid mixtures. For increased dipole moments m* A = m* B = 1.5 the corresponding 
schematic phase diagram is shown in Fig. d c )- ^ n this case Fig. EJshows the projections of the triple lines and critical 
lines onto a plane p = const. The vapor-liquid, vapor - ferromagnetic fluid, isotropic fluid - ferromagnetic fluid, and 
the A-rich ferromagnetic liquid - B-rich ferromagnetic liquid phase boundaries for dipole moments m* A = m B = 1.5 
are shown in Figs. El an dEl A further increase of the dipole moments to m A = m* B = 2 leads to the schematic 
phase diagram shown in Fig. d^)- A cut at Ap* = const through those coexistence surfaces is shown in Fig. 
describing the isotropic fluid ferromagnetic fluid and the A-rich ferromagnetic liquid B-rich ferromagnetic liquid 
phase transitions. 
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(4) We have also constructed the schematic phase diagrams of binary fluid mixtures of non-polar and various polar 
particles of equal size (see Fig. |S} and have discussed the topological changes of these phase diagrams due to an increase 
of the dipole moment of the polar component B. For dipole moments m\ — and m* B = 1 Fig. Elshows a cut through 
the schematic phase diagram (Fig. Efb)) describing the vapor- liquid, the A-rich isotropic liquid - B-rich isotropic 
liquid, and the vapor - ferromagnetic liquid phase coexistences. Figure 1101 displays another cut at A[i* = const 
through the phase diagram in Fig. Elb) indicating the possibility of A-rich isotropic liquid - B-rich ferromagnetic 
liquid coexistence besides the vapor-liquid one. Upon increasing the dipole moment of the polar component B from 
m B = 1 to m B = 1.5 the topology of the corresponding phase diagram shown in Fig. OUb) changes to that shown 
in Fig. [5Jc) which includes different types of phase coexistences: B-rich isotropic fluid - B-rich ferromagnetic fluid 
coexistence (see Fig. 1111 . A-rich isotropic fluid - B-rich ferromagnetic fluid coexistence (see Fig. I12p . and vapor - 
A-rich isotropic liquid coexistence together with A-rich isotropic liquid - B-rich ferromagnetic liquid coexistence (see 
Fig. I13|) . A further increase of the dipole moment to m* B = 2 leads to the phase diagram shown in Fig. |HJd) . 

(5) We have studied the dependence of phase diagrams of binary Stockmayer liquid mixtures on the size ratio 
obbIoaa- We have found that for 1 < ubb/^aa < 2 1 / 3 the phase diagrams of binary Lennard- Jones liquid mixtures 
(m A = m* B = 0) do not exhibit any topological change with increasing size ratio. For m A = m* B = Figs. 1141 and 
1151 demonstrate that with increasing particle size ratio only the thermodynamical parameter ranges of vapor- liquid 
and liquid-liquid coexistence are changed. In Fig. 1161 the schematic phase diagrams of binary Stockmayer liquid 
mixtures with particle size ratios o~bb/o~aa = 1 and o~bb/o~aa = 2 1 / 3 and for the dipole moment ratio m B /m* A = 2 
are shown. Such binary mixtures with dipole moment and particle size ratios different from 1 can be considered as 
a first approximation of polydisperse ferrofluids. The topological change of the phase diagram of binary Stockmayer 
liquid mixtures with dipole moments m* A = 0.75 and m* B = 1.5 due to the increase of the size ratio from <jbb/gaa = 1 
to obbI^aa = 2 1 / 3 is displayed in Figs. ll"oT a) and fToT c). Corresponding slices through these phase diagrams for 
Au* = const have been calculated (see Figs. 1171 and HS]l. Figure [T7I shows vapor - A-rich isotropic liquid and A-rich 
isotropic liquid - B-rich ferromagnetic liquid coexistence, while Fig. 1181 displays vapor - B-rich isotropic liquid and 
B-rich isotropic liquid - A-rich isotropic liquid coexistence. For dipole moments m\ = 1 and m B = 2 Figs. I16t b) and 
lltif d) show another type of topological change in the phase diagrams of binary Stockmayer liquid mixtures due to 
the increase of the size ratio from o~bb/o~aa = 1 to obbIo~aa = 2 1 / 3 . The corresponding slices through these phase 
diagrams for A/i* = const have been calculated (see Figs. 1191 and |2"U)) . Figure |2"U1 shows that with increasing size ratio 
vapor-liquid coexistence emerges besides the coexistence of B-rich isotropic fluid and B-rich ferromagnetic fluid (Fig. 
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